Surface-peaked effective mass in the nuclear energy density functional and its 

influence on single-particle spectra 
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Calculations for infinite nuclear matter with realistic nucleon-nucleon interactions suggest that 
the isoscalar effective mass of a nucleon at the saturation density, m* /m, equals 0.8 ± 0.1. This 
result is at variance with empirical data on the level density in finite nuclei, which are consistent 
with m*/m ~ 1. Ma and Wambach suggested that these two contradicting results may be reconciled 
within a single theoretical framework by assuming a radial-dependent effective mass, peaked at the 
nuclear surface. The aim of this exploratory work is to investigate this idea within the density 
functional theory by using a Skyrme-type local functional enriched with new terms, r(Vp) 2 and 
T^r, where r and p denote the kinetic and particle densities, respectively. We show that each 
of these terms can give rise to a surface peak in the effective mass, but of a limited height. We 
investigate the influence of the radial profile of the effective mass on the spin-orbit splittings and 
centroids. In particular, we demonstrate that the term quenches the I/5/2 — I/7/2 splitting in 
40 Ca, which is strongly overestimated within conventional Skyrme parametrizations. 

PACS numbers: 21.60. Jz,21.60.-n 
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I. INTRODUCTION 

Irrespective of their level of sophistication, infinite nu- 
clear matter (INM) calculations provide a relatively low 
isoscalar effective mass (IEM), m* /m « 0.8 ± 0.1, at the 
saturation density (TMHl • This finding is at variance with 
empirical data on the level density at the Fermi energy 
in finite nuclei, which are consistent with m* jm « 1. In- 
deed, the level density close to the Fermi surface appears 
to be very well described within simple models using bare 
nucleonic mass and a phenomenological one-body poten- 
tial, U(r), of, e.g., the Woods-Saxon form 

It is not at all obvious, however, whether the suc- 
cess of such non-self-consistent approach is built upon 
solid theoretical foundations or is a mere effect of the fit- 
ting strategy of U(r), which is adjusted precisely to the 
single-particle (SP) spectra. Indeed, rigorous and in prin- 
ciple exact treatment of the single-particle propagator 
(Green's function) is obtained within the Landau quasi- 
particle theory by solving the Dyson equation, which con- 
tains, in general, a spatially non-local and energy (fre- 
quency) dependent self-energy [ll|, E3- Hence, from the- 
ory point of view, the use of static, momentum and en- 
ergy (frequency) independent mean potentials seems to 
have no justification at all. 

Quantitatively, the most important effect on the SP 
spectra is expected to come from spacial non-locality or 
momentum dependence of the mean potential. This ef- 
fect is naturally taken into account in the Hartree-Fock 
(HF) approximation based on effective nucleon-nucleon 
(NN) interactions, like the most popular Skyrme [l3| 
or Gogny [l4[ forces. Indeed, the zero-range Skyrme 
(pseudo-) potential depends explicitly on momentum. In 
case of the finite-range Gogny interaction, the effect en- 
ters through the exchange term. The energy dependence 
is considered only as a perturbation correction atop of 



the HF result [15Hl9j and is usually completely ignored. 

The self-consistent HF method involving the Skyrme 
or Gogny interactions leads to a non-uniform, position- 
dependent IEM, m*(r)/m. Due to a rather restricted 
form of these forces, a typical IEM profile changes 
smoothly from the bare-mass value of 1 outside the nu- 
cleus to the INM value of 0.8 ±0.1 in the nuclear inte- 
rior. This lowers the nuclear level density at the Fermi 
surface usually well below the experimental value and, in 
turn, spoils the spectroscopic properties of these forces, 
see [2(3, HH and references quoted therein. 

The contradicting results concerning the SP level den- 
sity in INM and finite nuclei are related to the fact that 
in finite nuclei the SP levels couple to low-lying collective 
surface vibrations, which leads to an increase of the IEM 
at the surface. Several authors argued that this effect can 
be taken into account within the random phase approxi- 
mation, i.e., by going beyond mean field, see (l7l I22I [23jj . 
This viewpoint is difficult to reconcile with the density 
functional theory (DFT) that should warrant a proper ra- 
dial dependence of the IEM through the fit to empirical 
data provided that the energy density functional (EDF), 
or the effective NN interaction, is rich enough to accom- 
modate the anticipated increase of the IEM in the surface 
area. 

Taking into account the particle-vibration coupling 
may help resolve the discussed controversy in the follow- 
ing way. The IEM is smaller than one inside the nucleus, 
as required for nuclear matter. On the other hand, it is 
peaked up at the surface, so that it may give one when 
averaged over the nuclear volume. This may allow to si- 
multaneously reproduce the level density in finite nuclei. 

The surface-peaked IEM was first explored by Ma and 
Wambach. In their seminal paper (24[, they used a phe- 
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nomenological one-body Hamiltonian, 



For even-even nuclei, 



H(r) 



2m* (r) 



V + U(r), 



(1) 



consisting of the Woods-Saxon potential with the spin- 
orbit (SO) term, U(r), and the kinetic term with a radius- 
dependent IEM taken as follows. 



m 



dr 



(2) 



and depending on two adjustable constants, a and /3, 
respectively. The function g(r) is related to unperturbed 
self-energy, 



g(r) = S (r, k F (r), e F )/S o (0, k F (0), e F ). 



(3) 



It mimics the radial dependence of the SP potential and 
was chosen so that So( r , kp(r), ep) — U(r). It essentially 
follows the profile of the nucleonic density, g(r) ~ p(r). 
Ma and Wambach demonstrated that one can indeed 
unify the description of the SP levels in INM and finite 
nuclei by taking m*(r)/m 0.7 in the bulk and peaking 
it up at the nuclear surface to m*(r)/m ss 1.2 -j- 1.5, de- 
pending on the nucleus. This result was later confirmed 
by Farine et al. [25[ to hold in self-consistent HF calcu- 
lations. 

The aim of the present work is to explore the IEM mod- 
ifications caused by the spacial non-locality of the NN in- 
teraction and by the particle- vibration coupling within a 
fully self-consistent model. Unlike Farine et al. [25| , who 
obtained a surface-peaked IEM by adding new terms to 
the Skyrme interaction, we approach the problem from 
the perspective of the DFT, taking as a starting point the 
local Skyrme functional and adding new terms thereto. 
Nowadays, the nuclear EDFs are treated as basic enti- 
ties, independent of the NN interaction. In this sense, 
our approach is more general and more flexible than the 
Skyrme HF theory. 

In this exploratory work, we extend the Skyrme func- 
tional by adding two terms, which are described in de- 
tail in Sect. [TTJ Their effect on the radial IEM profile 
is presented in Sect. IIIII In Sect. IIV1 we examine their 
influence on the SP levels. Sect. [V] contains a discus- 
sion of the results, particularly in connection with the 
recent studies of the tensor term. The paper is concluded 
in Sect. I VII All considerations presented in this work 
concern the isoscalar effective mass, and are limited to 
doubly-magic N = Z nuclei. Such nuclei are spherical, 
which further simplifies the calculations. 



II. EXTENDED LOCAL ENERGY DENSITY 
FUNCTIONAL 

The local Skyrme EDF, TL(r), consists of the kinetic 
and potential parts, 



t=o,i 



(4) 



H t 



Po]Pt 



= ci 



Jt<2 



(5) 



YJ 



PtV ■ J t 



Index t = 0, 1 denotes isospin, and C t are adjustable 
coupling constants. The constant Cf[p ] is density- 
dependent and equals Cf [po] = Cf' + Cf' '(^p-) Q , where 
po,e 9 is the INM saturation density. The potential energy 
terms are bilinear forms of time-even particle, kinetic, 
and tensor densities, pt, rt, Jt, and their derivatives. The 
density 3 t is the vector part of the spin-current tensor, 
Jt,A = J2fj,v £\nv$t,nv Since we consider N — Z nuclei, 
only the isoscalar (t = 0) channel is active. Readers in- 
terested in details are referred, e.g., to Ref. [2a |. 

Within this formalism, the IEM comes from variation 
of H over to, 



^ _6H _ h 2 



(6) 



Thus, m*(r)/m depends on the particle density. It is 
equal to one outside the nucleus and, for positive values 
of Cq , smaller than one in the bulk. Therefore, the con- 
ventional Skyrme EDF contains only the volume term 
~ g(r) of parametrization ©. 

In order to allow for a surface-peaked profile of the 
IEM one has to add new terms to the functional. Having 
an almost complete liberty in choosing their form, we 
decided to add either the term 



^ A) (r)=^ (V " )2 r (r)(Vp (r)) 2 



or 



ui B \r) 



Crdpf dr 




ro(r) 



dpo(r) 
dr 



(7) 



(8) 



We call these two variants of the model A and B, and 
explore them separately, never including both terms si- 
multaneously. 

The rationale behind these two choices is almost self- 
evident. Both terms are proportional to the kinetic den- 
sity, To, and to derivatives of the particle density, po- 
Hence, their variations give 



2m* (r) 2m 



+ C^p (r) + C T (Vp) \\7 Po (r)f 



and 



2m* (r) 



— + C oPo (r) + C 



(9) 



(10) 



for variants A and B, respectively. They modify the 
Skyrme value ((6]) at the nuclear surface, where the deriva- 
tive of the particle density is large, and give small contri- 
butions elsewhere. Since the density gradient is negative 
in the surface region, while its square is positive, the peak 

is obtained for negative values of Cg*' Vp '' in variant A, 
and for positive values of (j^ d P/ dr m variant B. 
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The new terms also contribute to the mean potential 
according to the formula 



U (r) = 



8H 
Spa ' 



(11) 



This contribution is included self-consistently in our cal- 
culations. 

The term A is a scalar, so its form is natural even if 
spherical symmetry is not assumed. The use of term B 
is restricted to nuclei of spherical shapes, what seems to 
be a disadvantage. We have decided to explore this term 
mainly because it mimics closely bulk and surface terms 
w 1 - ag(r) + /3^M of the Ma and Wambach 
parametrization ([2]). 

The new terms introduced here are of completely dif- 
ferent nature than those used by Farine et al. [25| . In 
their model, the IEM is modified not by the gradient of 
the density, but through its value, 



2m* (r) 



2m 



CoVo(r) 



Cfr£ +1 (r), 



(12) 



where f3 is a new parameter. In such an approach, the 
peak at the surface is not due to the slope of the density, 
but because the density there equals approximately half 
of its value in the interior. 
In a preliminary study 



271 ] . we also investigated two 



other terms, r and rAp. These, however, do not pro- 
vide the expected radial profile of the IEM, and are not 
considered here. 



III. RADIAL DEPENDENCE OF THE 
ISOSCALAR EFFECTIVE MASS 

We base all our calculations on the SkXc Skyrme pa- 
rametrization of Ref . [28[ . This force was designed with 
particular attention to describe the SP levels, and has the 
IEM close to unity right from the beginning. We add our 
new terms to the SkXc functional, and vary some of its 
coupling constants as described below. We perform the 
calculations for three doubly-magic, N = Z, spherical 
nuclei, 40 Ca, 56 Ni, and 100 Sn. 

As already explained, we aim at obtaining an IEM pro- 
file which provides m*(r)/m lower than one inside the 
nucleus and peaked at the surface, so that its mean value 
is close to one. In order to fulfill the latter condition we 
require that 



rf J r 



poO) m*(r) 



.4 



m 



= 1, 



(13) 



where A is the mass number. A similar constraint was 
used by Ma and Wambach. Note that while the new 



coupling constants, Cq ' Vp ^ and C^ dp ^ dr , are responsible 
for the peak, Cq accounts for lowering the IEM inside the 
nucleus. Thus, the condition fp~3|) can be satisfied by an 
appropriate balance between the new coupling constants 
and Cq. 



TABLE I: Limiting values of the new coupling constants cor- 
responding to the largest peak and anti-peak in 40 Ca. 
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FIG. 1: Isoscalar effective mass, m* (r) /m, versus r in 40 Ca for 
variants A (left) and B (right). The solid (dashed) lines show 
the IEM profiles with the largest peak (anti-peak), obtained 
for the limiting values of the new coupling constants shown 
in Table [J 



We proceed with the calculations for each of the con- 
sidered nuclei separately in the following way. In both 
variants, we scan a range of the appropriate new cou- 
pling constant, and, for each value thereof, we readjust 
Cq to fulfill the constraint (fTB")) , keeping all the remaining 
coupling constants frozen at their SkXc values. 

For completeness, we examine both positive and neg- 
ative values of each new coupling constant, which corre- 
sponds both to a peak and to an anti-peak at the surface. 
It turns out that the range of each new constant is limited 
in both directions. Beyond that range, rapid oscillations 
in the density occur, and the HF iterations diverge. The 
limiting values of the coupling constants give rise to the 
largest peak or anti-peak that could be obtained. As an 
example, we give those limiting values for 40 Ca in Table 
HI and show the corresponding IEM profiles in Fig. [1] 

It can be seen from Fig. [T] that both variants of our 
model are rich enough to render a surface-peaked IEM. 
However, the height of the peak is rather limited. In 
variant A, the IEM in the nuclear interior is larger than 
^0.9 and reaches only ~1.10 at the surface. These values 
equal ~0.8 and ^1.25 for variant B, respectively. For the 
same nucleus, Ma and Wambach quote the values of ~0.7 
and ~1.25. This result can be reproduced only in variant 
B of our model. It should also be mentioned that we 
failed to obtain a surface-peaked IEM with those Skyrme 
parametrizations which have a small effective mass in the 
INM limit, like - 0.7 4- 0.8. 
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IV. INFLUENCE OF THE ISOSCALAR 
EFFECTIVE MASS GEOMETRY ON 
SINGLE-PARTICLE LEVELS 

We are now in a position to examine the influence of 
the new terms on SP spectra. Since our later discussion 
will mostly concern the SO properties, we consider the 
SO splittings, 



E so (n,£) = e{j = £-\,n,£)-e{j = £+\,n,£), (14) 
and the centroids of the SO partners, 



E c (n, £) = - [e(j = i - ±, n,£) + e(j = i + \, n,£)] 



. (15) 

Here, e(J,n,£) are the SP energies. Since we are inter- 
ested in delineating trends rather than details concern- 
ing the influence of the new terms on SP spectra we take 
e(J, n, £) simply as eigen-energies of the one-body Hamil- 
tonian, see discussion in Ref . 29] . To focus attention, we 
discuss only the neutron splittings and centroids. 

Let us remind that within the DFT formalism, the SO 
splittings are determined by the SO field, Bo(r), which 
is a derivative of T-L over Jo, and, for spherical symmetry, 
equals 



B (r) = C{j -Cj^, 



(16) 



where Jo stands for radial component of the vector part of 
the spin-current Jo- Since in most Skyrme parametriza- 
tions, including SkXc, the tensor coupling constant, C$ J , 
is very small or zero, the SO field is in practice propor- 
tional to the gradient of the density, which is largest at 
the surface. 

We perform two different types of calculations. In both 
cases we scan the whole admissible range of the new cou- 
pling constants, but: 

(i) For each value of the new coupling constant we 
readjust Cq to fit the constraint (|T5j) , as in Sect. 
IIII1 Hereafter, we call these calculations perturba- 
tive. 

(ii) In addition to adjusting Cq to the condition ([TU]) . 
we attempt to refit also other coupling constants 
to masses and radii of the considered nuclei. These 
calculations will be dubbed non-perturbative, and 
are described in details below. 



A. Perturbative calculations 

The SO splittings and centroids obtained from the per- 
turbative calculations in variant A are plotted in Fig. [5] 



in function of C, 



Note that the peak appears for 



negative values of this coupling constant. As already 
discussed, the term A modifies the IEM profile rather 
weakly. Consequently, it also affects the SO splittings 
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FIG. 2: (Color online) Spin-orbit splittings (upper panels) 
and centroids (lower panels) in function of C IT| - Vp - ) obtained 
from variant A of the perturbative calculations in 40 Ca (left), 
56 Ni (middle), and 100 Sn (right). 



and centroids in a relatively modest way. Exceptions are 
the I/5/2 — I/7/2 splittings in 40 Ca and 56 Ni, and the 
n = 2 centroids, 2s, 2p, and 2d. 

In order to understand these results, let us inspect 
Fig. |3l which shows the particle density, po, and its gra- 
dient in function of r for 40 Ca and 56 Ni, calculated in 

variant A for Cg equal to and to the limiting value 
of -280MeVfm 7 . It turns out that even the largest pos- 
sible peak only weakly affects the density profile in the 
surface region, but quite strongly amplifies the density 
fluctuations inside the nucleus. Thus, the gradient of the 
density and the SO field are only moderately increased 
at the surface, which does not influence the SO splittings 
significantly. On the other hand, the modifications in the 
interior have more impact on low-^ states with wave func- 
tions inside the nucleus, than on higher-^ states, whose 
wave functions are more spread toward the surface. 

The SO splittings and centroids obtained in variant B 
are presented in Fig. [U Note that here the peak appears 
for positive values of (j^ dp / dl \ The tendencies are very 
different than in variant A. With the development of a 
surface-peaked IEM, one observes a systematic and sub- 
stantial reduction of the SO splittings for high-^ orbits. 
In particular, the splittings for n = 1 orbitals, like Id, 
If, lg, change by as much as 2MeV throughout the en- 
tire range of Cg dp ^ dr . The splittings of the states with 
n = 2 and low— £, including 2s, 2p, 2d, are relatively 
weekly affected by the new term with no prevailing trend. 
However, these splittings are small, so their relative vari- 
ations can still be as large as 20%. The centroids of all 
the calculated orbits go steadily up with (j^ dp / dr \ There 
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is also a tendency to increase the distance, 



AE= \Ec(nJ)-E c (n + l,l-2)\ 



(17) 



between the SO centroids of, 
pseudo-spin partners. 



loosely speaking, the 



The changes of the SO splittings and centroids with 
(jT<ip/dr can ^ n i ce iy correlated with the changes in the 
radial dependence of the particle density shown in Fig. [5] 
With increasing (J rdp ^ dr , the density drops in the nuclear 
interior and the surface of the nucleus moves outward. 
Hence, also the gradient of the density and the SO field 
are repelled from the bulk, where the wave functions are 
located, and the SO splittings decrease. This argument 
is stronger for high— £, n = 1 orbits with an appreciable 
amount of the wave function in the surface region, and 
much weaker for low— £, n — 2 orbits, whose wave func- 
tions have one node inside the nucleus. The centroids 
move up in energy because the potential well becomes 
shallower and wider. 
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FIG. 4: (Color online) Spin-orbit splittings (upper panels) 
and centroids (lower panels) in function of (j Td p/ dr obtained 
from variant B of the perturbative calculations in 40 Ca (left), 
56 Ni (middle) and 100 Sn (right). 
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FIG. 3: Isoscalar particle density (upper panel) and its gradi- 
ent (lower panel) versus radius obtained from variant A of the 
perturbative calculations in 40 Ca (solid lines) and 56 Ni (dashed 
lines) for C ,T( - Vp - ) = (black lines) and for the highest-peak 
value of C T(Vp)2 = -280MeVfm 7 (gray lines). 



B. Non-perturbative calculations 

In the perturbative calculations, the new terms deteri- 
orate the performance of the SkXc functional concerning 
the INM saturation energy, binding energy per nucleon, 
and the incompressibility modulus, as well as radii, bind- 
ing energies, and SP levels in finite nuclei. The influence 
on the INM properties is illustrated in Fig. O This is so 
because the perturbative calculations take some coupling 
constants too far away from the SkXc values without 
readjusting the remaining ones. Hence, although helpful 
in understanding the response of the nucleus against the 
new terms, these calculations are of limited applicabil- 
ity. In this Section, we attempt to remove this drawback 
by refitting some other Skyrme coupling constants to se- 
lected radii, masses, and SP levels in finite nuclei. 

We freeze the tensor and SO coupling constants, Cq 

and Cq J , at their SkXc values. At each value of Cg^ Vp ' ) 
or CQ dp ^ dr , the constant Cq is adjusted to fulfill the con- 
dition (|13p . as in the perturbative calculations. Simulta- 
neously, Cq' , Cq ", and Cq Ap , are refitted to minimize 
the merit function 

x 2 = E - ~J l ] ( 18 ) 

i 1 

Here, wf h ^ and w^ xp%> denote the theoretical and exper- 
imental values of the z-th observable, and the weight of 
is an a priori assumed error associated with the i-th ob- 
servable. It may be understood as a desired theoretical 
accuracy. Observables used in the fit include masses and 
radii of 16 O, 40 Ca and 56 Ni, and the mass of 100 Sn. We 
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FIG. 5: Isoscalar particle density (upper panel) and its gradi- 
ent (lower panel) versus radius obtained from variant B of the 
perturbative calculations in 40 Ca (solid lines) and 56 Ni (dashed 
lines) for (j Td p/ dr = q (black lines) and for the highest-peak 
value of C rdp/dr = 70MeVfm 6 (gray lines). 



take of equal to 2MeV for masses and O.lfm for radii. 

It is not our intention to create a new functional here. 
We only aim at investigating the impact of the new terms 
on the SP levels while keeping the basic nuclear proper- 
ties close to the SkXc values. This is why we fix the 
coupling constants Cq and C^ J . Indeed, they influence 
the SO splittings and centroids directly through Eq. (IT51) . 
By fixing them, we can track the response of the split- 
tings and centroids solely to the new terms. 

In the non-perturbative calculations, the masses of 
40 Ca, 56 Ni and 100 Sn vary by less than ±1% through- 
out the whole range of the coupling constant (j^ dp ^ dr j n 
variant B. The mass of 16 varies within ±5%. The INM 
properties are shown in Fig. [5] For most values of the new 
coupling constants, the saturation density stays within a 
range of ±0.005/m" 3 (±3%) with respect to the SkXc 
value, the binding energy per particle changes by no more 
than ±0.6MeV (±4%) and the incompressibility modu- 
lus by less than ±5MeV (± 2%). These results guarantee 
that our refitted functional retains the basic characteris- 
tics of the SkXc parametrization and that new effects, if 
any, can be ascribed mostly to the new terms. The IEM 
profiles obtained from the non-perturbative calculations 
do not differ significantly from those given in Fig. [T] 

The neutron SO splittings and centroids are shown in 
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FIG. 6: (Color online) Saturation density (top), binding en- 
ergy per nucleon (middle), and incompressibility modulus 
(bottom) of infinite nuclear matter obtained from perturba- 
tive (green dashed lines) and non-perturbative (red solid lines) 
calculations in variant A (right) and B (left) in function of the 
corresponding new coupling constants. 
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FIG. 7: (Color online) Spin-orbit splittings (upper panels) 

and centroids (lower panels) in function of C T ( Vp ) obtained 
from variant A of the non-perturbative calculations in 40 Ca 
(left), 56 Ni (middle), and 100 Sn (right). 



Figs. [7] and [8] for variants A and B, respectively. In vari- 
ant A, there is almost no dependence of the splittings on 

C^ VpS> , and only a small increase in AE in 56 Ni and 
100 Sn occurring with the development of the peak. Vari- 
ant B is more promising. Indeed, the previously observed 
trend of splittings decreasing with the onset of the peak 
is retained. The centroids move slightly up in energy and 
the values of AE are becoming systematically smaller. 
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FIG. 8: (Color online) Spin-orbit splittings (upper panels) 
and centroids (lower panels) in function of C rdp ' dr obtained 
from variant B of the non-perturbative calculations in 40 Ca 
(left), 56 Ni (middle) and 100 Sn (right). 



V. DISCUSSION 



spin- unsaturated. 

It was demonstrated in Refs. [29|, |3l| that in or- 
der to reproduce the experimental value of the neutron 
I/5/2 ~~ I/7/2 splitting in 40 Ca, one has to dramatically 
reduce the SO strength, C t v,/ . This, in turn, requires 
a substantial increase of the tensor strength, C/ jf, in 
order to keep the SO splittings in spin-unsaturated nu- 
clei at a reasonable level. On the other hand, global fits 
to nuclear masses [32j do not allow for such significant 
modifications of the SO and tensor coupling constants. 

From this point of view, the reduction of the neutron 
I/5/2 — 1/7/2 splitting in 40 Ca by the surface-peaked IEM 
opens a possibility of reproducing this quantity with- 
out such a radical quenching of the SO strength. Con- 
sequently, also the tensor term would not need to be 
strengthened so much. This observations rise hopes that, 
with the inclusion of the surface-peaked IEM, the SP 
properties and nuclear masses may be both correctly de- 
scribed by a single functional. 

Another interesting observation is that also the split- 
tings between the pseudo-spin partners decrease with the 
emerging peak in variant B of the non-perturbative calcu- 
lations. These values are also overestimated by conven- 
tional Skyrme functionals, so that their reduction may 
improve the agreement with experiment. 



An important conclusion from our study concerns the 
behavior of the neutron I/5/2 — I/7/2 splitting in 40 Ca 
with the development of the surface peak in the IEM. 
This quantity is systematically overestimated in the con- 
ventional Skyrme functionals, independent of the para- 
metrization, see Ref. [29| and references therein. One can 

see, e.g., from Fig. @] for Cg dp ^ dr =0 that the SkXc func- 
tional misses the experimental value of 5.240MeV [3(| 
by more than 2MeV. From this point of view, the results 
obtained in variant A of our model are rather disappoint- 
ing because the concerned splitting increases even more 
with the onset of the peak in the perturbative calcula- 
tions, which is definitely undesired, and remains almost 
constant in the non-perturbative calculations. Variant B 
is by far more promising. It can be seen from Fig. @] that 
for (7q dr =70MeVfm 6 , the difference between theory 
and experiment drops below lMeV in the perturbative 
calculations. This decreasing trend is retained also in 
the non-perturbative calculations. 

These findings should be discussed in connection with 
the recent investigations of the tensor term in the Skyrme 
functional, see, e.g., Ref. 29, 31, 32]. The tensor term, 
C/Jt , together with the SO term, C^ J p t V-3 t , determine 
the SO field according to Eq. (|16|l and, consequently, the 
SO splittings. As explained, e.g., in Ref. [32|], the tensor 
term is inactive in spin-saturated nuclei, i.e., when the 
SO partners are both empty or filled, and is strongest 
in spin-unsaturated nuclei, i.e., when one SO partner is 
filled and one empty. From among the nuclei considered 
here, 40 Ca is spin-saturated, while 56 Ni and 100 Sn are 



VI. SUMMARY AND CONCLUSIONS 

We considered the SkXc Skyrme functional enriched 
with two new terms, C^ Vp ^ r (Vpo) 2 and Cq dp / dr TQ ; 
which we call A and B, respectively. They are designed 
to produce a surface-peaked radial profile of the isoscalar 
effective mass (IEM). Such a profile accounts for coupling 
of single-particle (SP) states to low-lying surface vibra- 
tions, and it may allow for a correct description of SP 
level density close to the Fermi surface both in infinite 
nuclear matter (INM) and in finite nuclei. 

It was demonstrated that both terms do indeed pro- 
duce a peak in the IEM, but of a limited height, because 
too large values of the new coupling constants lead to an 
instability and divergence of the HF iterations. 

It turned out that the response of the SP levels to the 
presence of the peak strongly depends on the variant of 
the model. While the SP properties hardly change in 
variant A, most of the SO splittings systematically de- 
crease with the development of the peak in variant B. 
This includes the neutron I/5/2 — I/7/2 splitting in 40 Ca, 
which is strongly overestimated by conventional Skyrme 
functionals. Thus, the inclusion of the new terms may 
help reproduce the empirical SO splittings. It may also 
be helpful in determining the correct strength of the ten- 
sor term, which is also partially responsible for the SO 
properties. Finally, the onset of the peak in variant B re- 
duces the energy distance between the pseudo-spin part- 
ners in 40 Ca, which also improves the agreement between 
theory and experiment. 
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Rigorous fits to bulk and SP data are necessary to 
corroborate or falsify our findings. However, this ex- 
ploratory work demonstrates that construction of lo- 
cal functionals by adding specific, physically well mo- 
tivated terms may constitute a promising alternative 
to t he g radient-expansion method proposed recently in 
Ref. (33J. 
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